Surface-induced heating of cold polar molecules 
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We study the rotational and vibrational heating of diatomic molecules placed near a surface at 
finite temperature on the basis of macroscopic quantum electrodynamics. The internal molecular 
evolution is governed by transition rates that depend on both temperature and position. Analytical 
and numerical methods are used to investigate the heating of several relevant molecules near various 
surfaces. We determine the critical distances at which the surface itself becomes the dominant source 
of heating and we investigate the transition between the long-range and short-range behaviour of 
the heating rates. A simple formula is presented that can be used to estimate the surface-induced 
heating rates of other molecules of interest. We also consider how the heating depends on the 
thickness and composition of the surface. 

PACS numbers: 34.35.+a 33.80.-b, 37.10.Mn 42.50.Nn, 



I. INTRODUCTION 

A number of techniques have recently been developed 
to cool polar molecules to low temperatures and to trap 
them for a second or longer. Using the switched electric 
field gradients of a Stark decelerator [l[ , polar molecules 
formed in a supersonic expansion have been decelerated 
to rest and then stored in electrostatic, magnetic or elec- 
trodynamic traps 0, H, 0, [i| . An electrostatic trap has 
been continuously loaded by filtering out the slowest frac- 
tion of the molecules present in an effusive beam Q . Po- 
lar molecules have also been cooled in a buffer gas of 
cold helium and then confined in a magnetic trap Q- 
Extremely cold polar molecules such as RbCs can be 
produced by the photoassociation of two species of ul- 
tracold atoms, followed by laser-stimulated state transfer 
Q. In all cases, the resulting molecules are typically far 
colder than their environments, and they may be heated 
by the absorption of blackbody radiation from that en- 
vironment. Unlike atoms, the polar molecules can be 
rotationally or vibrationally excited by their interaction 
with this blackbody radiation, and in many cases this 
can severely limit the trapping lifetime of the molecules. 
Indeed, the blackbody heating rate for trapped OH and 
OD has already been measured experimentally and found 
to limit the trapping lifetime to just a few seconds when 
the environment is at room temperature 0] . Calculations 
of the free-space heating rates for several polar molecules 
have already been presented [Toj . 

In most experiments so far, the cold polar molecules 
have been confined in macroscopic traps, with trap sur- 
faces typically several mm from the molecules themselves. 
There is now a great deal of interest in confining and 
manipulating these molecules much closer to surfaces, 
so as to build a 'molecule chip' technology analogous to 
that for atoms [ll| . Fast-moving molecules have already 
been trapped in travelling potential wells formed approx- 
imately 25 /im above a microstructured surface [l2j • This 
same structure can be used to decelerate the molecules to 
rest so that they can be trapped above the surface of the 



chip. Chip-based microtraps have been designed, along 
with schemes to interface the molecules with solid-state 
devices which could be used to cool, detect and control 
them coherently [HI]. Strong coupling to a supercon- 
ducting stripline cavity is possible when the molecules 
are just a few microns from the surface, and then the 
molecules can be the long-lived quantum memory of a hy- 
brid quantum information processor [T^ |. Recent devel- 
opments also herald the prospect of integrated molecule 
detectors, based on optical microcavities [151 ] or ultrathin 
optical fibers [l6[. 

These advances raise the question of the heating rates 
in the close vicinity of a surface. The influence of such a 
non-trivial environment on the internal atomic dynamics 
is commonly known as the Purcell effect [ItJ • Early theo- 
retical studies were devoted to the zero temperature case 
where the evolution is governed by spontaneous decay. 
As shown by linear response theory, the decay rate for 
an arbitrary environment can be given in terms of the 
classical Green tensor for the respective geometry [l8j . 
Alternative approaches have been developed on the ba- 
sis of classical electrodynamics pjl Ho[ and microscopic 
models [2l[ and have been applied to the case of an atom 
near a single surface or between two surfaces. Results for 
an arbitrary environment of electric [22j and magneto- 
electric bodies 23] , including local-field effects [24], HH , 
have also been obtained on the basis of macroscopic quan- 
tum electrodynamics (QED) and have been used to study 
atoms in bulk material [2^], outside (26[ or inside a mi- 
crosphere [24j . inside a spherical cavity 27] and even in 
the presence of left-handed meta-materials [U [2|| [29| . 
The linear-respons e ap proach has been generalised to fi- 
nite temperatures [301 ] where the internal dynamics is no 
longer governed by spontaneous decay alone, but stim- 
ulated emission and absorption of thermal photons also 
contribute. The respective environment-dependent tran- 
sition rates can again be expressed in terms of the classi- 
cal Green tensor; in addition, the thermal photon number 
comes into play. Ground-state heating rates of spinless 
atoms have been predicted to be very small near surfaces 
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[3l| . in contrast to the case of atoms with spin which 
have been investigated for planar surfaces, [3l|, HH wires 
[IH and carbon nanotubes [34| . 

In this paper, we calculate heating rates for a num- 
ber of polar molecules currently favoured by experi- 
menters. On the basis of macroscopic QED (presented 
in Sec. [TT|, we solve the internal molecular dynamics 
to obtain transition rates of a molecule in an arbitrary 
uniform-temperature environment (Sec. IIIip . In Sec. II VI 
the results are first used to calculate the rates in free 
space, and then as a function of distance from the sur- 
face of some common metals and dielectrics, as well as 
some unusual meta-materials. 



II. MACROSCOPIC QUANTUM 
ELECTRODYNAMICS AT FINITE 
TEMPERATURE 

Consider a molecule (or an atom) that is placed within 
an arbitrary environment of magneto-electric bodies. 
The coupled dynamics of the molecule and the body- 
assisted electromagnetic field can be described by the 
Hamiltonian (23l . l3a ] 



H = Ha + Hf 



H 



AF, 



where 



H A = J2 E n\n){r 



(1) 



(2) 



(E n , molecular eigenenergies; \n), molecular eigenstates) 
is the Hamiltonian of the molecule, 



J2 / d 3 r / dLonwfl(r,u)-f x (r,u), (3) 



H,- 



is the Hamiltonian of the electromagnetic field (including 
the internal charges present in the bodies) expressed in 
terms of the bosonic variables 



= = 



fl(r,u>)jUr',u>>) 



(4) 

Mr,w),fl tj (r',u/)\ =6 xx ,S ij 8(r-r')6(u;-u;'), (5) 



(note that f e is associated with the polarisation of the 
bodies and f m is related to their magnetisation) and 



Haf = d mn -E(r A )A r , 



(6) 



(d mn — (m\d\n), electric-dipole transition matrix ele- 
ments of the molecule; r A , molecular centre-of-mass po- 
sition; A mn — \m){n\, molecular flip operators) is the 
molecule-field interaction Hamiltonian in electric-dipole 



approximation. The electric field can be expressed in 
terms of the bosonic variables according to 



E{r) = \ du)E(r,u) +H.c, 



E(r,u)= J2 j dVG A (r,r».A(r», (8) 

\=e,m 

with the coefficients Ga being related to the classical 
Green tensor, G, by 



G e (r,r',u) = i 



2V -Ime(r',w)G(r,r',w), (9) 

C £ V 7T£q 



G m r, r', u = i - J | \' V'x G(r , r, lo)] . 

(10) 

For a given environment of macroscopic bodies, described 
by their linear, local and isotropic relative permittiv- 
ity e(r,(jj) and permeability p,(r, ui), the Green tensor 
is uniquely defined by the differential equation 



V x — V x - — e{r, u>) 

p(r,uj) c z 



G(r, r', u) = S(i — r') 

together with the boundary condition 

G(r,r',u)->0 for \r - r'\ oo. (12) 
The above definitions imply the useful integral relation 



£ J d 3 sG x (r,s^)-G+(r\s,u;) 



= ^ w 2 ImG(r,r',^). (13) 

7T 

In thermal equilibrium at uniform temperature T, the 
electromagnetic field may be described by the density 
matrix 



Pt 



c -H F /(k B T) 
tr [ c -iJ F /(fe B T)] 



(14) 



(fee, Boltzmann constant). Thermal averages (...) = 
tr[. . . fix] of the bosonic variables are thus given by 

(Mr,u>)) = = (ft(r,w)), (15) 

(/A(r,w)/v(r',<y)) = = </l(r, W )/X(r' )£l /)>, (16) 

(fl(r, u)f x ,(r', lo')) = n(w)6 xy S(r - r')6(u - J), 

(17) 

</ A (r,c,)/ A ,(r>')) 

= [n(w) + l]5 xx >6(r - r')5(u> - J) (18) 

where 



n(uj) 



J2 k ke~ khu) ^ kBl "> 



^ e -fcfiw/(fe B T) e huj/(k B T) _ i 



(19) 
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is the average thermal photon number. Recalling defini- 
tions (J7J) and |5]), the statistical properties of the electric 
field are found to be given by 

(E(r,uj)) = = (&(r,u)), (20) 
(E(r,u,)E(r',u')) = = (&{r,u;)&(r',u/)), (21) 



n(uj)u 2 ImG{r,r',oj)6(uj-Lu'), (22) 



[n(w) + l]w 2 Im G(r, r', - to') (23) 



where we have made use of the integral relation (p~3|) . 
Note that these relations are in accordance with the 
fluctuation-dissipation theorem [36j , 

( \ [M{r, u)E} (r', to') + E} (r', u')E(r, w)] ) 

= ^£ [ n ( w ) + i] u? Im G(r, r', w)<J(w - to'), (24) 

where the thermal photon energy is given by 



\hu> for k^T <C 
fce? 1 for k^T 3> 



^ [n(w) + 1] -> <; ;~ ~ 7; (25) 



in the zero- and high-temperature limits, respectively. 



III. INTERNAL MOLECULAR DYNAMICS 

Consider a molecule which is prepared at initial time 
t = in an arbitrary internal state, represented by its 
internal density matrix <r(0). The environment of the 
molecule is initially taken to be at uniform temperature 
T, so that the electromagnetic field is in a thermal state 
p(0)=p T . 

The internal molecular dynamics can be determined 
by solving the coupled equations 



i r°° 



(d nk A mk - d km A kn ) -E(rA, to) 



EJ (ta, to) • (d nk A mk - d km A k 



(26) 



and 



h(r,u) = -[H,f x {r,u)] 



-iojf\(r,Lo) + ^^2d mn -Gl(r A ,r,cu)A mn , (27) 



as implied by the Hamiltonian |T]) together with Eqs. ((2]), 
([3]) and ©. The electromagnetic field can be eliminated 



by formally solving Eq. ([27)) and substituting the result 
into Eq. (|26|) . For weak molecule-field coupling, the 
Markov approximation may then be employed to show 
that the dynamics of the internal density matrix of the 
molecule a is given by the equations (App. [X| 



0Vm(t) — — T n (7 nn (t) + ^ T kn (T kk (t), 



(28) 



&mn(t) — [— iw mB — q (r m + r n )] CT mn (t) 

for m 7^ n (29) 

(&mn = (m\a\n) = (A nm )). Here, the total loss rate T n of 
a level n is given by 



r„ = T n (r A ) = r n fc, 



(30) 



and the individual intra-molecular transition rates T nk 
from level n to level k read 

r„fc = T nk (r A ) = r° fe + r„ k 



2/Uo ~ 



^ u} nk d nk -luYG{r A ,r A ,\u) nk \)-d kn 
x{Q(u) nk )[n(u) n k) + 1] + 6(cj fc „)n(cD fcn )}(31) 
[0(z), unit step function] where 

r nfc = ^^Cul lk Q{u nk )d nk -lm^G(r A ,r A ,u nk )-d kn (32) 
and 

r «fe = ul k d nk -lmG(rA,r A , \ui nk \)-d kn 

x [Q(Q nk )n(Q nk ) + O(uj kn )n(co kn )] (33) 

denote the zero-point and thermal contributions to these 
rates [recall Eq. fl9"j)]. 

The intra-molecular transition rates depend on the 
shifted molecular transition frequencies 

Wmn = U mn (r A ) = LO mn + 5to m - 5u n , (34) 

where the frequency shift 

5u n = Suj n (r A ) = ^2 &Vnk, (35) 
fe 

of a given level n has contributions 
5u> nk = 6uj nk {r A ) = 5to° k + 5ioJ lk 

= ^ V J Q dej lu 2 ^d nk -ImG {1) (r A ,r A ,uj)-d kn 
n(oj) + 1 n{ui) 



LU nk - LO LO nk + LO 

Lo\d nk \ 2 \ n(uj) 



6irc 



n(uj) 



&nk — tO LO nk + LO 



(36) 
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(V, principal value) due to all other levels k, which can 
again be separated into their zero-point and thermal 
parts, 



<W° fc = Sco nk (r A ) 



2 d nk -ImG (1 \r A ,r A ,Lo)-d k 



dto LO 



L0 nk - LO 



(37) 



and 



&4fc = SuJ nk( r A) 

f-QG 

= -\V \ dujuj 2 d nk -lmG{r A ,r A ,Lu)-d k 

7™ Jo 



n(uj) 



+ - 



n(uj) 



^nk — LO LO nk + LO 



(38) 



respectively. Here, G*- 1 -* denotes the scattering part of 
the Green tensor according to the decomposition 



G(r,r',u) = G ia) (r,r',uj) + G (L) (r,r',Lo) 



(39) 



where the imaginary part of the bulk (free-space) part is 
given by '37] 



ImG (0) (r,r,Lo) 



6nc 



I 



(40) 



(1 5 unit tensor). The free-space zero-point frequency 
shifts associated with G^ ', i.e., the free-space Lamb 
shifts, are included in the bare transition frequencies io mn 
since they are determined experimentally in free space. 
The Green tensor being analytic in the upper half of 
the complex frequency plane, one can employ contour- 
integral techniques to rewrite the frequency-shift contri- 
butions as 



&L0 nk = Lo nk d nk -Re G {1) (r A ,r A ,{o nk )-d kn 



x {Q(6j nk )[n(Lo nk ) + 1] - <d(Lo kn )n(uj kn )} 
i 2/j, k B T 

N=0 



h 2 



(1 — $No)£,NOjjcn 



d n k-G^\r A , r A ^ N )-d kr , 



6ir 2 ch j 



du> lo 3 



N 

ti{lo) 



+ - 



n(u>) 



LOnk — LO LO nk + LO 



(41) 



[note that Re G(r, r' , —lo) = Re G(r, r', lo) for real to] with 
Matsubara frequencies 



2Trk B T 



N 



N, 



TV = 0,1,. 



(42) 



When neglecting the frequency shifts, the transition 
rates (|3TH) (f3"3"| obviously reduce to the well-known re- 
sults given, e.g., in Ref. [301 ] . 



It is worth noting that the internal molecular dynamics 
described by Eqs. (|2"8|) and ([2"9")) obeys probability conser- 
vation, 



tr Hi) = y^Annjt) 



^r nfc cr n „(t) + ^r kn <7 kk (t) = o, (43) 

n,k n,k 



dt 



where we have used Eq. (|30|) . From the property 

T nk = e h ^'^T kn (44) 

of the transition rates [see Eq. (f3"Tj) ]. it follows that in the 
long-time limit the molecule reaches a thermal state as 
its steady state 



&(t — > oo) = &t = 



tr [e 



J2 n E n \n){n\/(k B T) 
-J2 n E n \n)(n\/(k B T)] 



with 



E n = E n (r A ) = E n + HSlOn 



(45) 



(46) 



denoting the shifted molecular eigenenergies. This can be 
verified by noting that for this state the internal molec- 
ular evolution as given by Eqs. (|2"5|) and (|2"9"j) becomes 
static, 

&nn(t oo) = — r„tT„ n! T + ^ ^kn<^kk,T 

k 

= — ^nWnn,T 
k 

+ e- h ^ k ^ kBT ^r k n e h{D ™ k /( k ^a nn ^T = 0, (47) 

k 

a mn (t - oo) = e {-^-^ n]m{ t- t0)(TmntT = 
for m ^ n. (48) 

According to Eqs. (|30j) and (j3Tj) , the heating rate of 
a molecule prepared in its ground state |0) is given (ini- 
tially) by 



r = = 



T 
Ok 



2flQ 



^ H>l n(u>ko)d k-lm G(r A) r A) Lo k0 )-d ka , (49) 



due entirely to the absorption of thermal photons. 



IV. APPLICATIONS 

The energy associated with electronic excitation of 
molecules is typically large in comparison with thermal 
energy at room temperature, i.e. exp[— huo n Q / (k-BT)} -C 1, 
so according to Eq. (|4"4"|) . the fully thermalised state ef- 
fectively coincides with the electronic ground state. This 
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argument does not apply to the rotational and vibra- 
tional excitations of polar molecules, which occur at 
much lower frequencies. In this section, we study the 
ground-state heating rates Tq^ which provide a measure 
of the timescale on which this thermal excitation of the 
rotational and vibrational states takes place. We will 
assume that the frequency shifts induced by the envi- 
ronment are small enough to justify putting uj rnn =uj mn . 
In this case the thermal excitation rate from the ground 
state to state k becomes 



Ofc 



2^ 



ul a n(uj ko )d ok -ImG^fi, uj k0 ) -d k0 . (50) 



This has the great virtue that the temperature ap- 
pears only in the thermal photon number n(tu k o) [re- 
call Eq. (Unj)], while the position enters only through the 
Green tensor G. Therefore the dependence on temper- 
ature can be derived entirely from considering the free- 
space case, while the position-dependence can be under- 
stood completely from the behaviour at zero tempera- 
ture. 



electric dipole operator using Hund's case (a) basis states 
[38[. In this coupling scheme, the orbital angular momen- 
tum, L, is strongly coupled to the internuclear axis, and 
so is the electron spin, S, due to a strong spin-orbit cou- 
pling. The total angular momentum is J = L + S + R, 
where R is the angular momentum of the rotating nu- 
clei and is necessarily perpendicular to the internuclear 
axis. The projections of L, S and J onto the internu- 
clear axis are labelled by the quantum numbers A, E and 
fi = A + E. The projection of J onto the space-fixed z- 
axis is M. The basis states are labelled by the quantum 
numbers S, A, E, Q, J and M. 

For transitions between the rotational states, the ma- 
trix elements of the electric dipole operator are 

d mn = (QJM\d\n'J'M') = fi e {flJM\u\fl'j'M'} 

, ^ _j_ \ \ / \ -(- 1 \ 1 Q 

\ U mn ~ U mn)~f^ ' \ U mn ' U mn) /= ' U mn e z 



V2 



V2 



(55) 



where /i c is the molecular dipole moment at the equilib- 
rium internuclear separation, u = f/\r\, and 



A. Molecules in free space 

In free space, the Green tensor is given by Eq. (|4T))) . so 
the molecular transition rates become 

r _ r (0) _ r , pT 
J- nk — 1 n h — 1 nk > 1 nk 



\Unk\ 3 \d nk \ 2 

37rfeoc 3 



{Q(uj nk )[n(uj nk ) + 1] + 0(uj kn )n(uj kn )} 

(51) 



with 



and 



r 1 



F _ ^nfcl^fcl 2 Oil, , \ 



3|W . |2 



(52) 



[Q(LOnk)n(u>nk) + Q(uJ k n)n(uJ kn )] ■ 

(53) 

The total heating rate of a molecule initially prepared in 
its ground state thus reads 



3irhenc 3 



n(u) k0 ), 



(54) 



in agreement with Ref. [lOj . 

The ground-state heating rate of polar molecules will 
be dominated by transitions to the adjacent excited rota- 
tional and vibrational states, so we restrict our attention 
to these in the following. We calculate the heating rates 
for the set of ground state polar molecules listed in Tab.QJ 
which also gives the required molecular constants. 

We begin by considering rotational heating. To evalu- 
ate Eq. (|54[) we will calculate the matrix elements of the 



.9 

rnn 



= (-1) M -°V(2J 




(56) 



With this result, we obtain the selection rules for transi- 
tions between the basis states: Afi = 0, AJ = 0,±1, and 
AM = 0,±1. In this paper, we will not consider mixing 
of the electronic ground state with other electronic states, 
which leads to A-doubling, because the energy splitting 
that is induced is very small compared with the rotational 
energies and so does not alter any of our results. In this 
approximation, the states | ± f2 JM) are degenerate, and 
since Ail = we can confine our attention to the positive 
values of f2 only. While our equations make it clear how 
to handle initial states of given M', we will consider the 
initial molecular state to be unpolarised, averaging over 
the possible values of M'. 

The majority of the molecules listed in Tab. Q] have 
A = ground states. These molecules are best described 
using Hund's coupling case (b) [38[. The spin is not cou- 
pled to the internuclear axis and neither E nor f2 is de- 
fined. The rotational eigenenergies are 



E N = hB e N(N+l), AT = 0,1, 



(57) 



where B c is the rotational constant and AT is the rota- 
tional quantum number, N = J — S. The expansion of 
the E eigenstates in the case (a) basis is 38] 



\S,N,J,M)= J2 (-1) J ~ S V2N + 1 



Q=-S 



* n-n o |aJ ' M> (58) 
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Species 


Ground state 


B c (GHz) 


We(THz) 


Mc(KT 30 Cm) 


/ic(10" 21 C) 


m(10" 27 : 


LiH 


X X E+ 


222 


[39] 


42.1 


[10] 


19.6 [40] 


60.5 


[10] 


1.46 


NH 


X 3 E" 


500 


m 


98.4 


[41J 


5.15 [42J 


** 




1.56 


OH*" 


x 2 n 


555 


[43] 


112 


[44] 


5.56 [45] 


17.9 


[46] 


1.57 


OD*" 


x 2 n 


300 


[47] 


81.6 


[47] 


5.51 [45] 


*» 




2.97 


CaF 


X 2 E+ 


10.5 


[48] 


18.4 


m 


10.2 [49J 


172 


[41] 


21.4 


BaF 


X 2 E+ 


6.30 


[10] 


14.1 


[10] 


11.7 [10] 


285 


no] 


27.7 


YbF 


X 2 E+ 


7.20 


[50] 


15.2 


[51] 


13.1 [50] 


195 


[52] 


28.4 


LiRb 


X X E+ 


6.60 


[53] 


5.55 


[10] 


13.5 [54] 


21.4 


[10, 58] 


10.8 


NaRb 


X 1 E+ 


2.03 


[53J 


3.21 


[10] 


11.7 [54] 


12.6 


[10] 


30.0 


KRb 


X J E+ 


1.15 


[55] 


2.26 


[10] 


0.667 [54] 


1.89 


[10] 


44.3 


LiCs 


X X E+ 


5.80 


[56] 


4.92 


[10] 


21.0 [54] 


28.4 


no] 


11.1 


NaCs 


X 1 E+ 


17.7 


[56J 


2.94 


[10] 


19.5 [54] 


21.4 


[10] 


32.5 


KCs 


X X E+ 


92.8 


[56] 


1.98 


[10] 


8.61 [54] 


6.93 


[10] 


50.0 


RbCs 


X X E+ 


0.498 


[57] 


1.48 


[10] 


7.97 [54] 


4.41 


[10] 


86.0 



TABLE I: Properties of various diatomic radicals: electronic ground state, rotation and vibration constants, dipole moment 
and its derivative at equilibrium bond length, and reduced mass. For comparison with the constants used in [l(J, see [58| ]. 
* Reduced masses are given on the basis of the atomic masses (most abundant isotopes) of the molecular constituents as stated 
in Ref. gl]. 

** For NH and OD, the electric-dipole matrix elements for the transition between ground and first excited vibrational states 
can be given as |d fc | = 1.80 x 10~ 31 Cm 42] and \d ok | = 7.54 x 10" 32 Cm [§(}, respectively. 

*** The spin-orbit coupling constants required for OH and OD are ,4 = -4.189 THz and A = -4.174THz [glj, respectively. 



Using Eqs. ([55)) , l[56")) and ([58)1 . summing over the possible 
final states and averaging over initial states of different 
M', we find J2 k \ d ok\ 2 = l4 for X S, 2 £, and 3 S molecules. 
For 2 £ molecules, the ground state |iV = 0, J =1/2) can 
be excited either to \N=1, J=l/2) or to \N = 1, J=3/2), 
with branching ratios 1/3 and 2/3 respectively. The spin- 
rotation interaction lifts the degeneracy between these 
states, but this splitting is very small and we do not 
need to include it. For 3 E molecules, the ground state 
\N = 0, J = 1) can be excited to the three states with 
N = 1 and J = 0, 1, 2, with branching ratios 1/9, 1/3 and 
5/9 respectively. Again, we can neglect the small spin- 
rotation interaction that lifts the degeneracy between the 
three states. 

The electronic ground states of OH and OD are 2 II 
states and, for low values of J, are best described using 
Hund's coupling case (a). The Hamiltonian describing 
the fine structure contains a rotational part and a spin- 
orbit coupling, H fs = hAL-S + hB c (J — L — S) 2 . The 
rotational term couples states of the same J but different 
|f2|. Writing the matrix elements of the Hamiltonian as 
wfi,n' = (flJM\H ts \n'JM) we have dU 

^3/2,3/2 = ±hA / 2 + hB c [J(J + 1) - 3/4 T 1], (59) 

1/2,1/2 

ma/2,1/2 = -hB oy /(J + 3/2)(J- 1/2). (60) 

Diagonalizing this Hamiltonian gives a pair of energy 
eigenvalues for each value of J > 1/2, 

Ej = hB e [(J + 1/2) 2 - 1 ± Q/2], (61) 

where 

Q = V 4(J + 1/2) 2 + A/B C (A/B C - 4) (62) 



We will use the labels Fi and F2 to denote the states 
of lower and higher energy, respectively. For the low- J 
levels of OH and OD, the mixing of fl states is small be- 
cause |A| is considerably larger than B e J. Recalling that 
A is negative for these molecules, we can then identify Fi 
as having predominantly 2 n 3 / 2 character, and F2 as pre- 
dominantly 2 n!/2. For J = 1/2 there is only one level, 
which is of pure fl — 1/2 character. The eigenstates are 



[-Pi, J, M) =c+(J)|l/2, J, M) + c_(J)|3/2, J, M), 

J = 3/2,5/2,..., (63) 
\F 2 J,M) =c+(J)|3/2, J,M) -c_(J)|l/2, J,M), 

J = 1/2,3/2... (64) 



where 



c±(J) = y/1/2 ± (A/B c - 2)/(2Q) . (65) 



Using the selection rules between the basis states, we 
see that the possible transitions out of the molecular 
ground state \F\,J = 3/2), are those to the states 
(a)|F 1; J = 5/2), (b)|F 2 , J=l/2), (c) \F 2 , J = 3/2) and 
(d) \F 2 ,J = 5/2). Applying Eqs. ([55]) and d56)) to each 
of these four transitions, summing over the M sublevels 
in the final state, and averaging over the M' sublevels in 
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the initial state, we obtain 

£ |do.| 2 = [§4(3/2)4(5/2) + |c 2 _(3/2)c 2 _(5/2) 

k(a) 



6 1 



+ 5V3 

k(b) 

Z2 \ d okf 

fc(c) 



c + (3/2) C+ (5/2) C _(3/2)c_(5/2) ^ (66) 



1 „2 
3 



4(3/2)M e 2 



A<(3/2)d(3/2)^, 



(67) 
(68) 



k(d) 



10fc| 



4(3/2)c 2 _(5/2) + §4(3/2)4(5/2) 
3 c + (3/2) C _(3/2) C+ (5/2)c_(5/2)l Ai 2 . (69) 



With these preparations, we can now evaluate the rates 
for free-space rotational heating out of the ground-state, 
for the molecules listed in Tab. [J The lifetimes, 
= (r(°))- 1 , arc given in Tab. [IT] for environmental tem- 
peratures of 293 K and 77 K. Since there is little vari- 
ation of the dipole moment, the lifetime is mainly de- 
termined by the power of the thermal spectrum at the 
transition frequency. Apart from the weakest transi- 
tions in OH and OD all these lines lie on the low side 
of the peak frequency in the thermal spectrum, which 
is 17THz at 293 K or 5THz at 77 K. Note that the ro- 
tational constant is roughly given by B c w H/(4:wmR^) 
where R e is the equilibrium internuclear separation and 
m = m,\m<il {yn\ + 772,2) is the reduced mass, so as a rule 
of thumb, rotational heating is most severe for the light 
molecules. Strong heating is seen for LiH, NH, OH, OD, 
whose lifetimes are in the range of 2-6 seconds. For KCs 
and NaCs the heating is much less severe, and for the 
rest it is negligible for most practical purposes. Tab. |TT] 
also shows that rotational heating of OH and OD is dom- 
inated by transition (a), with the other transitions pro- 
viding small corrections to the heating rate, even though 
they are at higher frequencies. This behaviour is due to 
the exceedingly small transition dipole moments of the 
latter transitions. The rotational excitation lifetimes of 
all these molecules can be extended by going to lower 
environmental temperatures. Figure [T] illustrates this 
temperature-dependence in the light molecules LiH, NH, 
OH, OD and KCs. 

Let us next turn our attention to vibrational heating. 
To a good approximation, the vibrational eigenenergies 
of the deeply-bound states of a diatomic molecule are 



huj e {v + |), v = 0, 1, 



(70) 



where uj c is the vibrational constant and v is the vibra- 
tional quantum number. The corresponding eigenstates 
are 



1/4 



1 



: H v (y/aq)e 



(71) 



Species 


^(GHz) 


1 j 2 
l^k \ d Ok\ 

«! 


r (0) 

293 K 


(s) 

77 K 


LiH 


444 


1 


2.1 


9.1 


NH 


999 


l 


6.4 


31 


OH 






2.1 


17 


(a) 


2.51 xlO 3 


0.405 


2.4 


18 


(b) 


3.80xl0 3 


0.00999 


49 


550 


(c) 


5.64xl0 3 


0.00775 


34 


720 


(d) 


8.67xl0 3 


0.00124 


120 


8,400 


OD 






6.3 


37 


(a) 


1.41 xlO 3 


0.402 


7.2 


39 


(b) 


3.93xl0 3 


0.00381 


120 


1,400 


(c) 


4.89 x10 s 


0.00302 


110 


1,800 


(d) 


6.48 x10 s 


0.000636 


340 


10,000 


CaF 


21.0 


1 


3,400 


13,000 


BaF 


12.6 


1 


7,200 


28,000 


YbF 


14.4 


1 


4,400 


17,000 


LiRb 


13.2 


1 


4,900 


19,000 


NaRb 


4.05 


1 


70,000 


260,000 


KRb 


2.30 


1 


6.7 xlO 7 


2.5 x 10 s 


LiCs 


11.6 


1 


2,600 


10,000 


NaCs 


35.5 


1 


330 


13,000 


KCs 


186 


1 


62 


250 


RbCs 


0.995 


1 


2.5 xlO 6 


9.5 xlO 6 



TABLE II: Free-space lifetimes for rotational heating out 
of the ground state at 293 K and 77 K. For OH and OD, the 
effects of the transitions (a)-(d) (see main text) are also shown 
separately. Also given are the frequency and the square of the 
dipole matrix element for each transition. For comparison 
with the results of LLOl, see 15811. 




T(K) 



150 200 300 



FIG. 1: Free-space lifetimes of the ground state against rota- 
tional heating as a function of environment temperature for 
LiH (thick solid line), NH (thick dashed line), OH(a) (thick 
dotted line), OD(a) (thin solid line) and KCs (thin dashed 
line). 



where q = R — R , R being the internuclear separation, 
H n are the Hermite polynomials and a — 2irm(v e /h. Ex- 
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v, t~ii"»/~' i no 
op L-UlCb 


2tt l 1Mz i 


_(0) 
T 


\ S ) 

rp — 771/- 


LiH 


42.1 


25 


6.5 xlO 9 


NH 


98.4 


310,000 


1.3 xlO 25 


OH 


112 


9.8 xlO 6 


2.2 xlO 29 


OD 


81.6 


200,000 


3.7 xlO 21 


CaF 


18.4 


4.7 


23,000 


BaF 


14.1 


1.8 


1,300 


YbF 


15.2 


4.1 


4,700 


LiRb 


5.55 


128 


2,700 


NaRb 


3.21 


1,400 


13,000 


KRb 


2.26 


120,000 


850,000 


LiCs 


4.92 


80 


1,300 


NaCs 


2.94 


580 


4,900 


KCs 


1.98 


12,000 


74,000 


RbCs 


1.48 


63,000 


350,000 



TABLE III: Lifetime against free-space vibrational heating 
out of the ground state for various polar molecules at 293 K 
and 77 K. For comparison with the results of [JtJ, see [58l ]. 

panding the electric-dipole operator in a Taylor series 
about the equilibrium separation, d w (/x c + fj,' e q)u, and 
recognizing that the first term cannot couple different 
vibrational states, we write the matrix elements for rovi- 
brational transitions in the form 

(vVtJM\d\v Vl'JM') = (i' e {QJM\u\Cl'j'M'){v\q\v'). 

We see from this equation that the rovibrational tran- 
sitions must satisfy the same rotational selection rules 
as already given above, and that to leading order in q, 
the vibrational selection rule is Av = ±1. For transitions 
between v' = and v = 1 we have 

We neglect the contribution of rotational energy to the 
transition frequency since it is typically smaller than the 
vibrational energy by two orders of magnitude. This 
means that we can simply add up the contributions of 
transitions (a)-(d) in calculating the transition dipole 
moments for OH and OD. Thus we obtain, 

k 

where / rot = 1 for the S molecules, while for molecules 
with a 2 n 3 / 2 ground state, 

/ rot = #4(3/2) + §c 2 _(3/2) + A c 2 +(3 /2)ci(3/2). (75) 

The calculated lifetimes for free-space vibrational heat- 
ing out of the ground state are given in Tab. [Till f° r 
T = 293 K and T = 77 K. These lifetimes are mainly de- 
termined by the vibrational transition frequencies. Since 



T<°)(i) 




100 150 200 250 300 

FIG. 2: Ground state vibrational heating lifetimes in free 
space vs temperature for LiH (thick solid line), CaF (thick 
dashed line), YbF (thin solid line), and BaF (thick dotted 
line) . 



uj c oc 1/ x/m, the lightest molecules have the highest vi- 
bration frequencies, which lie above the 17THz peak 
of the room temperature spectrum, whilst the heaviest 
molecules vibrate well below this frequency. The vibra- 
tional transition frequencies of CaF, BaF and YbF fall 
close to this maximum, and of the molecules considered 
these three also have the largest values of /i' e . For both 
reasons, the ground-state lifetimes of these molecules are 
limited by vibrational heating to less than 5 s. For LiH 
and LiCs the vibrational heating is an order of magni- 
tude slower, whilst it is exceedingly slow for all the other 
molecules. This slowness is mainly due to inefficient cou- 
pling with the thermal radiation which occurs both for 
the heavy molecules LiRb, NaRb, KRb, NaCs, KCs and 
RbCs whose vibration frequencies are too low and, even 
more strikingly, for the light molecules NH, OH and OD 
whose frequencies are too high. Due to the large transi- 
tion frequencies, the impact of lowering the environment 
temperature is even more striking for vibrational heating 
than in the rotational case. This is illustrated in Fig. [2] 
where the temperature-dependence of the lifetimes is dis- 
played for the molecules LiH, CaF, BaF and YbF which 
are most strongly affected by vibrational heating. 

The relative importance of rotational vs vibrational 
heating varies from molecule to molecule. Rotational 
heating dominates for the hydrides and for NaCs and 
KCs, while vibrational heating is dominant for the fluo- 
rides in the list, and for the other alkali dimers. 

We have confined our attention to the rates for rota- 
tional and vibrational excitation out of the ground state. 
The calculation is very easily adapted to the excited 
states, remembering that then there will be both exci- 
tation processes to higher lying states, and de-excitation 
processes to lower lying ones. For the latter processes, 
n(u)ko) should be replaced with n(ujka) + 1 to account for 
spontaneous emission. The calculations also need to be 
modified if applied electric or magnetic fields are present, 
so as to account for the Stark or Zeeman shifted transi- 
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tion frequencies, and any associated change in the tran- 
sition dipole moments. 



B. Molecules near a surface 

We turn now to the question of how proximity to a 
surface can affect the heating rate. Let us consider a 
molecule at distance z A from the surface of a homoge- 
neous magneto-electric body of (relative) permittivity 
s(lo) and (relative) permeability /J,(to). The body can be 
modelled by a semi-infinite half space provided it is close 
enough to the molecule and sufficiently smooth. The 
scattering part of the Green tensor is then given by 



G (1) (r,r,w) 



1 

8^ 



dq | e 2 ^ 



P 2 c 2 



p ) fex^-x H~ &y&y 



where 



mMI - Pi 



9^! 



gW - Pi 
e[u)P + Pi 



(76) 



(77) 



are the reflection coefficients for s- and p-polarised waves, 

(78) 



P = 



q 2 , 



Pi 



(Im P, Im Pi > 0) denote the z-component of the wave 
vector in free space (/3) and inside the half space (Pi) 
and q is its component parallel to the surface. For com- 
putational purposes, it is often convenient to express the 
Green tensor as an integral over /3, 



G (1) (r,r,w) 



1 
8^ 



1 

8^ 



dpe 2l>3z 
+ 2\1 
dbe~ 2bz 



P 2 c 2 



p ) fex^-x H~ &y&y 



P 2 c 2 



TpG z G z 



b 2 C 2 

2 J K^X&X ~T~ CyBy) 



2 1 



b 2 c 2 



(79) 



(P = ib). Here, the first term represents the oscillating 
contributions due to travelling waves, while the second 
term contains the exponentially decaying contributions 
from evanescent waves. 

Transition rates for a molecule near a half space can be 
obtained by substituting the scattering part of the Green 
tensor [as given by Eq. ([76]) or (175)) ] together with its 
free-space part [Eq. P0)) ] into Eq. (|31|) . The ground-state 



heating rates then take the particularly simple form 

r(r A ) = r(°)+r«(r A ) 



r (o) 



2ttc 

^kO 



Im trG (1) {r A ,r A , 



<-^ko) 



(80) 



In general, the integral appearing in Eq. ([76]) or ([75]) has 
to be evaluated numerically, but analytic results can be 
obtained for sufficiently small or large molecule-surface 
separations. The non-retarded limit applies to short dis- 
tances, where z\^/~Ep,\tjj / c <C 1, while the retarded limit 
holds for long distances such that zw/c> 1. For the ma- 
terials we consider in this paper, ly 7 ?/!! takes on values 
between 65 and 27,000 for rotational transitions and be- 
tween 4.5 and 700 for vibrational transitions, depending 
on the molecule and the material. Consequently, there is 
quite a large range of intermediate distances where nei- 
ther limit applies. 

In the non-retarded limit, the Green tensor ([79]) is 
dominated by the integral over evanescent waves, which 
effectively extends up to a wave vector 6= l/(2z). Over 
most of this region, P ~ Pi ~ iq, allowing us to use the 
approximations 



1 



e{u) + 1 



(81) 



Performing the remaining integral and retaining only the 
leading order in zlu / c, one finds that in this non-retarded 
limit, the Green tensor is well approximated by [3l| 



G (1) (r,r,w) 

c 2 e{u) - 1 



32ttw 2 z 3 s(uj) + 



\&%&x ^"y^-y ~t~ 2g 2 6z)' (82) 



Note that by retaining only the leading order in zuj/c, 
the dependence on r s and thus also that on /i has van- 
ished. In any case, fi is close to 1 even for the ferromag- 
netic metals at the typical frequencies of interest here (i.e. 
uj/2tt > 10 GHz). On substituting Eq. §2§ into Eq. (pITj) . 
we obtain the approximate, near-field transition rate 

r ( \ - r (0) -u \ dnk \ 2 + \ d »k.z\ 2 Ime(u nfc ) 

i nk [Z A) -l nk+ gTOo ^3 | £Kfe) + 1 |2 

x {@(Lu nk )[n(Lu nk ) + 1] - Q(uj kn )n(bj kn )}. (83) 

In particular, the ground-state heating rates ([80[) are 
approximated by 



T(z A ) = r(°) (i + f&) 



where 




2|e(w w ) + l| 



(84) 



(85) 
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is a scaling length that applies to calculations in the non- 
retarded limit. For a metal with permittivity 



Material ajp(rad/s) 7(rad/s) oj P /7(rad/s) 



1 



lu(lo + ry) 



(86) 



and provided that the transition frequency is sufficiently 
small, uJko ^ 7 < w p , z m may be estimated by the simple 
relation 



7 



(87) 



The plasma frequency, lup, and damping constant, 7, are 
given for various conductors in Tab. IIVI 

We stress that Eq. ([84]) applies only in the non- 
retarded limit, and that the distance z nr typically lies 
well outside this limit. We define a second relevant length 
scale, z c , the characteristic distance at which the surface- 
induced rate becomes equal to the free-space rate. This 
does not coincide with z nr , because the non-retarded limit 
is not valid at this distance. We have calculated z c by nu- 
merical integration of Eqs. ([?9"[) and ([8H)) and we present 
the results for molecules near a gold surface in Tab. [Vj 
and for a range of other conductors in App. [Bj The ta- 
ble also gives the corresponding values of z nr , which are 
typically 2-5 times smaller. 

Since Eq. ([8"4"[) does not apply at length scales in the 
vicinity of the critical distance, we searched for an al- 
ternative formula by fitting to the numerical results ob- 
tained from the integration of Eqs. (|79|) and ([80]) at dis- 
tances z < z c . We find that for the molecules and sur- 
face materials studied, the heating rates throughout this 
range are well approximated by the empirical formula 

r(*) = r<«»(i + 4 + 4)- 

V Z A Z A / 

Furthermore, a fit to the set of critical distances for rota- 
tional heating given in App.lBl suggests the approximate 
formula 



7 



3c 



4 V 2ul u\ Q 



(89) 



This empirical formula was found to be accurate to 
within 1% for all the surfaces and molecules considered, 
except in cases where the critical distances are particu- 
larly small (the hydrides and KCs), where deviations be- 
tween 1% and 10% are more typical. The same formula 
does not accurately predict the critical distances for vi- 
brational heating, but these are of less importance due 
to their very small values. We stress again that Eq. (f88|) 
is only empirical as the z7 2 term has no physical inter- 
pretation. 

The critical distances given in Tab. [V] show that the 
surface does not generate any significant heating when 
the molecules are more than a few hundred /jm away. 
However, if the molecules are held a few /im from a sur- 
face, as they might be on a molecule chip, there is a 



Au 


1.37 xlO 16 


4.12 xlO 13 


4.53 xlO 18 


Al 


2.25 xlO 16 


1.22xl0 14 


4.15 xlO 18 


Pd 


8.36 xlO 15 


2.16xl0 13 


3.24 xlO 18 


Ag 


5.77xl0 15 


1.15xl0 13 


2.89 xlO 18 


Cu 


1.12xl0 16 


4.41 xlO 13 


2.87xl0 18 


Mo 


1.14xl0 16 


7.86 xlO 13 


1.65 xlO 18 


Fe 


6.23xl0 15 


2.79xl0 13 


1.39 xlO 18 


Co 


1.18xl0 16 


1.07xl0 14 


1.29xl0 18 


W 


9.72 xlO 15 


8.53 xlO 13 


1.11 xlO 18 


Ni 


7.44 xlO 15 


6.53 xlO 13 


8.49 xlO 17 


Ft 


7.75 xlO 15 


1.04xl0 14 


5.75 xlO 17 


ITO 


3.33 xlO 15 


1.68 xlO 14 


6.63 xlO 16 



TABLE IV: Drude parameters for various conductors. Values 
are taken from Ref. [621. with the exception of those for ITO 
(indium tin oxide) [63|] ■ The list is in order of decreasing 
uJp/7, which corresponds to increasing surface heating rate. 



Rotational 



Vibrational 



Species 


Znr (/mi) 


z c (fim) 




Zc(H 


LiH 


0.73 


1.9 


0.035 


0.071 


NH 


0.42 


1.0 


0.020 


0.042 


OH 






0.018 


0.0039 


(a) 


0.23 


0.50 






(b) 


0.17 


0.36 






(c) 


0.13 


0.27 






(d) 


0.10 


0.20 






OD 






0.022 


0.0048 


(a) 


0.34 


0.78 






(b) 


0.17 


0.35 






(c) 


0.15 


0.30 






(d) 


0.12 


0.24 






CaF 


5.5 


19 


0.061 


0.12 


BaF 


7.8 


27 


0.072 


0.14 


YbF 


7.1 


25 


0.069 


0.14 


LiRb 


7.6 


26 


0.13 


0.27 


NaRb 


17 


64 


0.19 


0.41 


KRb 


24 


98 


0.24 


0.54 


LiCs 


8.2 


29 


0.15 


0.30 


NaCs 


3.9 


13 


0.21 


0.44 


KCs 


1.3 


3.6 


0.27 


0.60 


RbCs 


42 


180 


0.33 


0.75 



TABLE V: Non-retarded length scales and critical distances 
for surface enhancement of rotational and vibrational heating 
rates near a gold surface. 



substantial increase in the rotational heating for all the 
molecules considered, apart from the hydrides. Even in 
cases where the free-space rate is small, the enhanced 
rate can be very large because of the rapid inverse-power 
scaling. For example, in free space, the rotational heat- 
ing time of CaF, 3400 s, is enormous compared with the 
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1000 r 



(a) 




o.i L — 
1 

z c {jim) 
5.00 r 



5 10 50 100 5001000 



w fc0 /(2^)(GHz) 



(b) 




uW(2tt)(THz) 



FIG. 3: Exact critical distances for surface-induced heat- 
ing vs frequency of the molecular transition, (a) Rotational 
heating, (b) Vibrational heating. Surface materials are gold 
(circles), iron (squares), platinum (diamonds) and ITO (trian- 
gles) . Frequencies of the plotted data points correspond (left 
to right) to RbCs, KRb, NaRb, LiCs, BaF, LiRb, YbF, CaF, 
NaCs, KCs, LiH, NH, OD(a) and OH(a). Solid lines indicate 
the slope corresponding to the lj~ 3 ^ 4 frequency dependence 
given by the empirical formula (|89[l . 



vibrational lifetime of 4.7 s. However, at a distance of 
1 fun from a room temperature gold surface the lifetime 
for rotational excitation drops to about 8 s and at smaller 
distances the rotational heating rate dominates over the 
vibrational rate. For the hydrides, the high rotational 
frequency that gives them rapid free-space heating also 
makes them relatively insensitive to the proximity of the 
surface except at sub-micron distances. 

Figure [3](a) shows the critical distances for rotational 
heating of various molecules near a range of surfaces. It 
is seen that their frequency scaling follows quite nicely 
the a>~ 3 / 4 dependence given by Eq. (fB"9")) , which is indi- 
cated by the solid line. This trend continues in Fig.^b), 
which shows the critical distances for vibrational heating. 
These are, of course, smaller because the vibrational fre- 
quencies are higher. 

Equation (|87|) shows that the short-range heating de- 
pends on the surface material through the factor wp/7. 
A low value of this ratio leads to a large critical distance 
and hence to more surface-induced heating. The values 
are displayed in the last column of Tab. IIVI for various 
metals in order of decreasing Wp / 7. In this list, gold is the 
metal of choice when trying to minimise surface-induced 
heating of the molecules, as also indicated by the circles 



in Fig. [3] At the opposite extreme lies indium tin oxide 
(ITO), which we include here because of its attractive 
combination of conductivity and optical transparency. It 
has a low plasma frequency and a high damping rate 7 
and so generates stronger heating, as shown by the tri- 
angles in Fig. [3] The values for other metals generally 
obey the Wp/7 scaling, though there are some exceptions 
where values of Wp/7 are very similar. 

As indicated by Eq. ([55]) . the surface- induced heating 
will be particularly large in cases where |e| is not large, 
but e has a significant imaginary part. This never hap- 
pens for a conductor, but can occur for dielectric media 
that happen to be strong absorbers at the relevant fre- 
quency. As an example, consider borosilicate glass which 
has £ = 6.2 + 0.21i for frequencies in the tens of GHz range 
(64j. The critical distance for rotational heating of CaF 
near such a surface is about 620 fim, very much larger 
than for a metallic surface. The timescale for rotational 
heating, which is 3400 s in free space, is thus reduced to 
just 0.14 s when this molecule is 10 fim from such a glass 
surface. 

Next, we turn to the retarded limit, where zuj/c^> 1, so 
the integrand in Eq. ([75)1 is rapidly oscillating or decaying 
over most parts of the integration regime. The main 
contribution to the integral (|79p comes from the region 
around the stationary-phase point q — 0, so that we may 
approximate 



(90) 



The integral can then be performed, and upon retaining 
the leading order in c/(zw), one finds that the Green 
tensor in the retarded limit reads 



G {1 \r,r,uj) = 



8ttz v /^7) + v /i(w) 

x (e x e x + e y e y ). (91) 



Consequently, the transition rates (|3"Tj) are given by 

r I \ - r>(°) 1 W nk\\^nk,x? + \dnk,y\) 
1 nk{ZA) — 1 nk H - — 



ATT£Qhc 2 ZA 



x Im 



\/ H(c0nk) + \/e(u) n k) ) 

x {B(u} nk )[n(uj nk ) + 1] - 0(uj kn )n(uj kn )}; (92) 
for a good conductor, |e| 3> \fi\, they further simplify to 



(0) ^l k (\ d rLk,x\ 2 + \ d nk,y\ 2 ) . ( ' 2z A UJ nk 



T nk {z A )=T nk ~ 



sin 



^ An£ hc 2 ZA "*"V c 

x {Q(uj nk )[n(uj nk ) + 1] - Q{uj kn )n(ujkn)}- (93) 

In particular, the ground-state heating rates (|80[) are 
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given by 



T(z A ) = r(°) 



1 + 



2zAOJko 



xlm f VmKo) - y^Ko) ^iz^o/c 



r (o) 



1 



2zALOkO 



2zA(jJnk 



(94) 



Thus, the surface-induced modification of the heating 
rates in the retarded limit is an oscillating function of 
distance, where the amplitude of the oscillation follows a 
power law. In particular, the heating rates approach 
their free-space values in the limit za — > oo. 

In order to see the entire distance dependence of the 
heating rate it is necessary to calculate the rates as given 
by Eq. ([79^1 and ([50]) numerically. The results are dis- 
played in Fig. Q] where we show the total heating rates 
as a function of distance for OH, LiH, CaF and NaCs 
molecules at distances in the range 1 — 500 /an from a 
gold surface. In all cases, the vibrational heating rate is 
dotted, the rotational rate is dashed and the total rate is 
a solid line. For OH the heating is entirely dominated by 
the rotational transitions over the whole of this distance 
range, and the vibrational contribution does not even 
appear in the plot. The heating rate is modulated with 
a period of 60 /im just as expected in the retarded limit 
[Eq. I|94p]. The heating rate is not greatly altered from its 
free space value, even at the shortest distance considered. 
For LiH, the heating is again dominated by the rotational 
transitions. The far-field oscillations modulate the rate 
and we see roughly one cycle with a period of 338 /tm. 
The heating rate rises sharply inside the critical distance 
for rotational excitation, which is 1.9 /im, whereas the 
vibrational contribution, having a much shorter critical 
distance, remains essentially constant down to 1 /im. For 
CaF, the heating is dominated by vibrational excitation 
at 18.4 THz, corresponding to an oscillation period in the 
far field of 8 /im, which can clearly be seen. Inside the 
19 /im critical distance for rotational heating, we see a 
dramatic increase in the rotational contribution to the 
rate, such that the two contributions are roughly equal 
at a distance of 1 /im from the surface. For NaCs, the two 
contributions are roughly equal in the far field and both 
are rather small. The 3 THz vibrational heating exhibits 
the expected far-field oscillations, whilst the rotational 
heating is at too low a frequency to show oscillations 
over this range. Inside the 13/im critical distance, the ro- 
tational heating increases rapidly, becoming a thousand 
times faster at a distance of 1 /im. 

In Fig. [51 we show once again the heating rate for NaCs 
as a function of distance from a gold surface (solid line). 
This figure also shows for comparison the heating rates 
near iron and ITO surfaces. At distances large enough 
for the retarded limit to apply, the heating rate given by 
Eq. ([M| is independent of the particular metallic surface 
since these are all good conductors at the relevant ex- 
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FIG. 4: Heating rates for OH, LiH, CaF and NaCs vs distance 
from a gold surface. Solid lines: total heating rate. Dotted 
lines: vibrational excitation rate. Dashed lines: rotational 
excitation rate. 



citation frequencies. At short range however, where the 
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FIG. 5: Heating rate for ground state NaCs as a function 
of the distance from gold (solid line), iron (dashed line) and 
ITO (dotted line) surfaces. 



FIG. 6: Heating rate for ground state LiH as a function of 
distance from fictitious meta-materials with e(uko) = fj,(uko) 
= ±10+i (solid line), e(u>ko) = 10+i, /i( w fco) = — 10+i (dashed 
line) and e(u>ko) = —10+i, fi(uJko) = 10+i (dotted line). 



near-field limit of Eq. applies, the heating rate be- 
comes proportional to j/ujp- As shown in Tab. IIVI this 
ratio differs widely between these materials and is a hun- 
dred times larger for ITO than for gold. For this reason, 
the ITO surface produces a larger heating rate at short 
distance and exhibits a longer critical distance than gold, 
as seen in Fig. [5j 

So far, we have discussed surface-enhanced heating in 
the presence of metallic and dielectric surfaces. It is also 
interesting to consider the heating rate for molecules in 
the vicinity of meta-materials, since these offer tunable 
magneto-electric properties [65, 66], and can even be left- 
handed [13] • A left-handed medium is realized when the 
real parts of both e and \i are simultaneously negative, 
leading to a negative index of refraction and a number 
of counterintuitive effects associated with the propaga- 
tion of the electromagnetic field inside such a medium 
[63] ■ Since the surface-enhanced heating rate of a sin- 
gle interface depends solely on the reflected electromag- 
netic field, one would expect it to be insensitive to left- 
handedness. For weakly absorbing media, the oscilla- 
tions seen in the retarded limit are small when the signs 
of Ree and Re/z are both positive or both negative. In 
the case where Re e = Re /z, and the imaginary parts are 
small, there are no oscillations at all since r s and r p are 
then very close to zero. The amplitude of the oscilla- 
tions is greatest when these reflection coefficients have 
their maximum values of 1. As we have already seen, 
this occurs for metals because |e| is much larger than 
In the context of meta-materials, reflection coeffi- 
cients close to unity are obtained for any weakly absorb- 
ing medium where Ree and Re/i have opposite signs. 
Note that this result is insensitive to the magnitudes 
of Ree, Re/i, which neither need to be equal nor par- 
ticularly large; they need only be of opposite sign and 
considerably larger than the imaginary parts. A meta- 
material engineered with these properties would produce 
large oscillations in the heating rate with a phase deter- 
mined by the chosen values of Re e, Re fi. Figure [S] shows 
the heating rate of a LiH molecule near fictitious weakly 



absorbing meta-materials with Ree = ±Re/i. We see 
that the left-handed material (Ree, Re/i<0) gives rise 
to exactly the same heating rate as a comparable ordi- 
nary material with Ree, Re/z>0 (the two curves can- 
not be distinguished on the plot), and that the oscilla- 
tions are suppressed. On the other hand, media with 
Re e = —Re^ result in large oscillations of the heating 
rate in the long-distance regime, with a phase that de- 
pends on the material properties. 



C. Surfaces of finite thickness 

The results of the previous section have shown that 
metallic surfaces can considerably enhance surface-indu- 
ced heating. In the context of chips, metal surfaces are of- 
ten unavoidable since they are used in current- or charge- 
carrying structures. One possible strategy to reduce the 
associated molecular heating is to reduce the thickness 
of the metal substrates. For a slab of finite thickness 
d, coated onto an infinitely thick substrate of permittiv- 
ity e s = e s (w) and permeability fj, s = fi s (w), the surface- 
induced heating rate is still given by Eq. (|80f together 
with Eq. ([75]) . but the reflection coefficients are now given 
by 



H 2 [3f3 s - /Xg/3f + UMfh [Ms/3 - ft] cot(/?iri) 
/i 2 ftft + Msft 2 + [Ms/? + ft] cot (ft d) 
e 2 ftft - e s ft 2 + ieft [e s/ 3 - ft] cot(ftrf) 
e 2 ftft + e s ft 2 + ieft [e s/ 3 + ft] cat (ft d) 



where 



ft 



■ e s (uj)[i s (LL>) - q 2 



(95) 
(96) 

(97) 



(lmft>0). 

Let us first consider the influence of the metal surface 
alone by letting e s = fi s = 1. In the non-retarded limit, 
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FIG. 7: Heating rate of NaCs as a function of distance 
from an ITO surface with thickness 10 /im (thick solid line), 
1 /im (dashed line), 0.1 fim (dotted line) and 0.01 /im (dashed- 
dotted line). The thin solid line shows how the result changes 
for a 0.01 /im thick surface when the ITO is coated onto 
borosilicate glass. 



A reduction of the short-range heating rates below the 
values observed for thick plates is eventually found but 
only once the coating is unfeasibly thin. Thus our cal- 
culations show that the heating at short-range cannot 
be reduced by reducing the material thickness. To un- 
derstand the increase of the heating rate with decreasing 
thickness, note that in the non-retarded limit the heating 
is mainly due to the coupling of the molecule with the 
surface plasmons at the front face of the plate. As the 
thickness decreases, these couple to the plasmons at the 
back face of the plate, leading to mutual enhancement 
and thus to an increase of the heating rate [68|]. To in- 
clude the borosilicate glass substrate we took e s (uJko) = 
6. 2 + 0. 21i for the rotational transitions [64[ and e s (uJko) — 
6.4 + 0.74i for the vibrational ones (69|. We find identical 
results whether or not this substrate is included, except 
for the thinnest coating, d = 0.01 /im, where we find that 
the presence of the substrate slightly reduces the heating 
rate, as shown by the thin solid line in the figure. 



the reflection coefficients may then be approximated by 



SUMMARY AND CONCLUSIONS 



H 2 + 1 + 2/i coth(gd) 



e 2 -l 



e 2 + 1 + 2e coth(gd) 



(98) 
(99) 



[recall the discussion above Eq. (fBTj) ]. Since q< 1/(2za), 
the short-range heating rate will be identical to that of an 
infinitely thick plate provided d^> za, since in this limit 
the above reflection coefficients reduce to those given in 
Eqs. ([77)) . On the other hand, the reflection coefficients, 
and hence also the heating rate, must become very small 
when d\e\ <§C za- We note immediately that, for molecule- 
surface separations of interest, a conducting surface needs 
to be unfeasibly thin for this limit to be reached, because 
of the enormously large values of |e| for a conductor. The 
behaviour between the two limits has to be determined 
from a numerical analysis. 

In the retarded limit, one may approximate 



e + (i + 2i y/ejl cot (^/ejlud/c) 



(100) 



so for a good conductor, |e| ^> \fi\, the reflection coeffi- 
cients and the heating rate become independent of the 
plate thickness at long range. 

In Fig.[jj we display the surface- induced heating rate of 
a NaCs molecule near ITO plates of various thicknesses as 
a function of molecule-plate separation. At large separa- 
tions, the heating rates are independent of the thickness, 
as predicted from Eq. (|100p . Over the entire range of dis- 
tances calculated, the ITO plate of thickness 10/mi gives 
the same result as a plate of infinite thickness. Reducing 
the thickness below this value increases the heating rates 
at short distances, contrary to the expectation of reduced 
rates at short range anticipated from Eqs. (|9"8")l and (|99p. 



Using macroscopic QED, we have calculated the inter- 
nal dynamics of a molecule placed within an arbitrary 
environment of magneto-electric bodies of uniform tem- 
perature. The internal time evolution is governed by 
the molecular transition frequencies and transition rates 
which depend on position and temperature. The depen- 
dence on temperature is due to the thermal photon num- 
ber and can be completely understood from considering 
the free-space case, while the position-dependence, which 
enters via the classical Green tensor for the electromag- 
netic field in the environment, can be derived from the 
behaviour at zero temperature. 

We have used the general formulae to study the 
ground-state heating rates of several polar molecules of 
current experimental interest, as a function of the dis- 
tance from various surfaces. We have given a simple ap- 
proximate formula that can be used to estimate the heat- 
ing rates for any other molecules at any distance from 
any surface of interest. For light molecules, particularly 
the hydrides, rotational heating dominates and limits the 
free-space lifetime of the ground state to a few seconds 
when the environment is at room temperature. For the 
metal fluorides we studied, vibrational heating dominates 
and again the room temperature free-space lifetime of the 
ground state is of the order of a few seconds. When the 
molecules approach a metallic surface, the heating rate 
can be greatly enhanced. This is particularly true for the 
rotational transitions where the critical distance at which 
the surface dominates the free-space rate is typically in 
the 1 — 100 /im range. For the hydrides, the free-space 
heating rate is large because the rotational frequencies 
are large, but this same fact also means that the criti- 
cal distance for surface-induced heating is rather small. 
Therefore, these molecules could be trapped up to a few 
microns from a surface with little change in the heating 
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rate. The effect of the surface on rotational heating is 
very much stronger for the heavier molecules, but since 
the rate in free-space is typically very small for these 
molecules, they too have lifetimes of a second or more at 
distances up to 1 /*m from the surface. 

We have shown that, at long range, the heating rates 
become insensitive to the particular surface properties, 
while at short-range the heating is faster for smaller val- 
ues of the parameter Wp/7. Of the metals considered, 
gold induces the lowest heating rate. We have also shown 
that decreasing the thickness of the surface tends to in- 
crease the heating rate, particularly at short distances. 
Dielectric materials that are good absorbers at the rele- 
vant frequency result in large critical distances and hence 
very large heating rates at short range. 

In the context of molecule chips, where confinement 
of molecules a few microns from the chip surface is en- 
visaged, our work shows that surface-induced heating 
should be considered carefully when selecting appropriate 
molecules and surfaces, but that confinement for several 
seconds in single quantum-states is quite feasible even 
when the surface is at room temperature. In all cases, 
lowering the environment temperature will allow for even 
longer lifetimes. For approach distances smaller than 
1 /jm, surface-induced heating becomes rapidly problem- 
atic, and cooling to cryogenic temperatures may be re- 
quired. 
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APPENDIX A: MARKOV APPROXIMATION 



Substituting the formal solution 



f x (r,LU,t) = C -^f x (r,u;) 



/ dre- iw ( t - r )d mn .G3;(r Aj r,a;)i tBn (r) (Al) 

n m * JO 



to Eq. (f27| into Eq. (|26| and using the integral rela- 
tion ([TBI, one obtains 



- » i f°° 

A mn (t) = iu mn A mn (t) + T X! / dW 



~ iuit [d nk A mk (t) - d km A kn (t)] -E{r A ,u) 

e lult E\r A ,uj)-[d nk A mk {t) - d km A kn {t)] 

+ Z„ m (t) (A2) 



with 



x d n k-lmG(r A ,rA,u)-dij 
- [e-^-^i^i^r) -e i ^ t -^A lj (r)A nk {t)] 

x d km -lmG(r A ,r A ,u;)-dtj} (A3) 

denoting the zero-point contribution to the internal 
molecular dynamics. This differential equation can be 
solved iteratively by substituting the self-consistent solu- 
tion 



A m nit) — e 1 """ 1 * A mn (0) 



illf "»/> 



Jo 



x je 1UT [d nk A mk (T) - d km A kn (T)]-E(r A ,Lu) 

+ e ilAjT E}{r A ,uj)-[d nk A mk (T) - d fem i fe „(r)] j (A4) 

to the truncated Eq. (|A2|1 without Z mn (t) back into 
Eq. (|A2|) . where at this level of approximation, the op- 
erator ordering in Eq. (|A4[) may be chosen arbitrarily. 
Taking expectation values according to Eqs. (p~5|) — (p~8|) . 
one arrives at 

iuJ mn (A mn (t))+(f mn (t))+(Z mn (t)) (A5) 

where 



(T mn {t)) 

= - £1 |J J™ dco Lu 2 n(iu) jf * dr [e-W-r) + 

x {e i ^ fc (*- r )[(i m; (r))d„ fc TmG(r A ,r yl ,^)- ( i H 

- (Aik (r) ) d nk Tm G(r A , r A , u) ■ d/ m ] 
- e^C^' [(i w (r))d km .ImG(r J 4,r A , w).d ni 

- (Ai n (T))d km -lmG(r A ,r A ,u)-di k ]^ (A6) 



16 



denotes the thermal contribution to the internal molec- 
ular dynamics. For weak molecule-held coupling, the 
contributions (T mn (t)') and (Z mn (t)') may be evaluated 
by means of the Markov approximation, i.e. by letting 
(j4 mn (r)) ~ e~ lu ' mn ( t ~ T '(A mn (t)') and evaluating the re- 
maining time integrals according to 

Jo 

~ it5{uj - Cb mn ) + VP - — . (A7) 

Assuming the system to be non-degenerate, so that off- 
diagonal molecular density matrix decouple from each 



other as well as from the diagonal ones, we arrive at 
Eqs. (HD and together with Eqs. (3DjH|38J). 



APPENDIX B: CRITICAL DISTANCES 

We have calculated the critical distances for surface- 
induced enhancement of ground-state heating rates for 
various combinations of molecules and materials on the 
basis of the data given in Tabs. Uland llVI The results are 
summarised in Tab. IVIl 
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.17 


.23 


.28 


.25 


.31 


.32 


.33 


.38 


.69 
























CaF 


19 


20 


21 


21 


24 


25 


26 


26 


28 


31 


53 


.094 


.074 


.071 


.086 


.11 


.098 


.13 


.13 


.13 


.15 


.30 


BaF 


28 


30 


31 


31 


35 


37 


37 


39 


42 


46 


78 


.12 


.093 


.088 


.11 


.14 


.12 


.16 


.16 


.16 


.20 


.37 


YbF 


25 


27 


28 


28 


32 


33 


34 


35 


38 


41 


71 


.11 


.087 


.083 


.10 


.13 


.12 


.15 


.15 


.15 


.18 


.35 


LiRb 


27 


29 


30 


30 


34 


35 


36 


38 


40 


44 


76 


.27 


.22 


.20 


.26 


.32 


.39 


.36 


.36 


.38 


.43 


.78 


NaRb 


66 


70 


72 


72 


82 


86 


88 


91 


97 


110 


180 


.42 


.37 


.34 


.43 


.51 


.48 


.55 


.57 


.60 


.68 


1.2 


KRb 


100 


110 


110 


110 


130 


130 


130 


140 


150 


160 


280 


.55 


.51 


.47 


.57 


.68 


.65 


.73 


.75 


.79 


.89 


1.6 


LiCs 


30 


32 


33 


33 


37 


39 


40 


41 


44 


49 


83 


.30 


.25 


.22 


.29 


.36 


.32 


.39 


.40 


.42 


.48 


.85 


NaCs 


13 


14 


14 


14 


16 


17 


17 


18 


19 


21 


36 


.45 


.40 


.37 


.46 


.55 


.52 


.59 


.61 


.64 


.73 


1.3 


KCs 


3.7 


3.9 


4.0 


4.0 


4.6 


4.8 


4.9 


5.1 


5.5 


6.0 


10 


.61 


.58 


.54 


.64 


.75 


.73 


.81 


.83 


.88 


.99 


1.7 


RbCs 


190 


200 


210 


210 


240 


250 


250 


260 


280 


310 


530 


.76 


.74 


.70 


.81 


.95 


.94 


1.0 


1.1 


1.1 


1.2 


2.2 



TABLE VI: The exact critical distances z c (fim) for the enhancements of ground-state heating rates of various polar molecules 
due to the presence of surfaces of various materials. 
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